The structural distortion in antiferromagnetic LaFeAsO investigated by a 

group-theoretical approach 
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As experimentally well established, undoped LaFeAsO is antiferromagnetic below 137K with the 
magnetic moments lying on the Fe sites. We determine the orthorhombic body-centered group 
Imma (74) as the space group of the experimentally observed magnetic structure in the undistorted 
lattice, i.e., in a lattice possessing no structural distortions in addition to the magnetostriction. 
We show that LaFeAsO possesses a partly filled "magnetic band" with Bloch functions that can 
be unitarily transformed into optimally localized Wannier functions adapted to the space group 
Imma. This finding is interpreted in the framework of a nonadiabatic extension of the Heisenberg 
model of magnetism, the nonadiabatic Heisenberg model. Within this model, however, the magnetic 
structure with the space group Imma is not stable but can be stabilized by a (slight) distortion of 
the crystal turning the space group Imma into the space group Pnn2 (34). This group-theoretical 
result is in accordance with the experimentally observed displacements of the Fe and O atoms in 
LaFeAsO as reported by Clarina de la Cruz et al. [nature 453, 899 (2008)]. 
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I. INTRODUCTION 

Undoped LaFeAsO undergoes an abrupt structural dis- 
tortion from tetragonal to orthorhombic [l|-|3| or to mono- 
clinic (^] symmetry at IbbK as well as an antiferromag- 
netic spin ordering transition at ^ [H-Il]- Clarina 
de la Cruz et al. \^ studied with high accuracy the struc- 
tural distortion of LaFeAsO at AK by neutron-scattering 
experiments and found the Fe and O atoms to be slightly 
shifted from their positions in the tetragonal phase. In 
agreement with the space group PA/nmm, the z coordi- 
nates of iron and oxygen are exactly z — 1/2 and z — 0, 
respectively, in the tetragonal phase. The values Az of 
the displacements at AK in ±z direction are reported to 
Az = 0.0006 and Az = 0.0057 for the iron and oxygen 
atoms, respectively These displacements cannot be 
simply explained by magnetostriction. 

In the framework of a nonadiabatic extension of the 
Heisenberg model [3, called nonadiabatic Heisenberg 
model (NHM) pia], the magnetic order and the low- 
temperature structural distortion in LaFeAsO may be 
understood by group-theoretical methods. While the 
non-magnetic distortion of LaFeAsO between \bbK and 
137K shall be discussed in a following paper in the con- 
text of the superconducting state in the doped material, 
the present paper investigates the magnetic structure and 
the associated structural distortions below 137K. 

In Sec.|llT]we shall determine the space group S of the 
experimentally observed magnetic structure in LaFeAsO 
under the assumption that there are no distortions in ad- 
dition to the magnetostriction. We shall show that S is 
compatible with the symmetry of the Bloch functions of 
the energy band of LaFeAsO denoted in Fig.[T]by the bold 
line. In the following Sec. llVl we will ask whether S may 
be the space group of a stable magnetic structure, and 
in Sec.|V]we shall present the concept of "allowed space 



groups" . Finally, in Sec. I VII we shall propose that the 
experimentally observed distortions of the crystal corre- 
spond to a change from the non-allowed space group S 
of the undistorted crystal to an allowed space group. 

The NHM does not distinguished between orbital and 
spin moments. Therefore, we always speak of "magnetic 
moments" which may consist of both orbital and spin 
moments. 



II. COORDINATE SYSTEMS 

The coordinates systems used in this paper are de- 
picted in Figs. [5] and [3] First, we use the x,y,z co- 
ordinate systems as depicted in Fig. [5] which coincides 
with the coordinates normally used in the literature on 
LaFeAsO, as, for instance, in the tables of Ref. [4]. This 
coordinate system defines the point group operations in 
the way described in Fig. [21 

In the following we shall consider three structures in 
pure LaFeAsO whose unit cells together with the basic 
translations Ti , T2 , and T3 of the corresponding Bravais 
lattices are given in Fig.[3l Fig. [3] (a) shows the tetragonal 
primitive Bravais lattice F^ of paramagnetic LaFeAsO 
above 155K which has the space group P4/nmm 
[15} . In Fig. [3] (b) the experimentally observed anti- 
ferromagnetic structure in a hypothetically undistorted 
crystal is depicted, which has the orthorhombic body- 
centered Bravais lattice rjj. Finally, Fig. [21(c) shows the 
antiferromagnetic structure in distorted LaFeAsO below 
137 K with the orthorhombic primitive Bravais lattice Fq. 

For the description of these three structures we need 
to resort to three coordinate systems with unusual rela- 
tionships. While the basic translations Ti, T2, and T3 
are adapted to the respective structure, the origin of the 
X, y, z coordinate systems and the directions of its axes 
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FIG. 1. Band structure of tetragonal LaFeAsO as calculated by the FHI-aims program [3,[T3, with symmetry labels determined 
by the authors. The symmetry labels can be identified from Table HI see Appendix. The bold line shows the magnetic band 
(as defined in Sec. IIII.2|l consisting of two branches. 



relative to the atoms are fixed for all the three structures 
in Figs. [3] Therefore, we may compare directly the space 
groups of the structures because the same point group 
operation in any structure is notated by the same sym- 
bol. As a consequence, the A and B axes (as given in 
Fig. but not the x and y axes, are perpendicular to 
each other in the orthorhombic structures. Furthermore, 
this choice of the x,y, z axes requires a renaming of the 
point group operations given in the tables in the text- 
book of Bradley and Cracknell [l^, as described in the 
notes to Tables H HI andHnl 



III. THE UNDISTORTED MAGNETIC 
STRUCTURE 



III.l. The experimentally determined magnetic 
structure and its magnetic group 



We show that the group Imma = ^Z^^ (74) with the 
orthorhombic body-centered Bravais lattice F^ depicted 
in Fig. [3] (b) is the space group of the experimentally de- 
termined magnetic structure as it is given by Fig. 4 of 
Ref. [4] . In Fig. [3] (b) the unit cell of antiferromagnetic 
LaFeAsO is depicted where only the Fe atoms are indi- 
cated since the magnetic moments are localized on the 
Fe sites [H, This Fig. H (b) is identical with Fig. 4 of 
Ref. [3| because by application of the basic translations 
Ti , T2 , T3 to the four Fe atoms in the unit cell the mag- 
netic structure in Fig. 4 of Ref. Q can be constructed. 



In this section we show that the experimentally ob- 
served magnetic structure in LaFeAsO is compatible with 
the symmetry of the Bloch functions of the magnetic 
band denoted in Fig.[T]by the bold line. First, in the fol- 
lowing subsection we shall determine the magnetic group 
of the experimentally observed magnetic structure under 
the assumption that there are no structural distortions 
in the material in addition to the magnetostriction. 



In a first step we show that the positions of the atoms 
of LaFeAsO are invariant under the symmetry operations 
of Imma. They are clearly invariant under the transla- 
tions of Imma since its basic translations Ti,T2,T3 as 
given in Fig. [3] (b) are also translations of the paramag- 
netic structure depicted in Fig. [3] (a) [cf. the following 
relation ^] . In order to show that the atoms of LaFeAsO 
are also invariant under the rotations and reflections of 
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FIG. 2. Cartesian coordinate system with the A and B axes 
defining the point group operations as used in this paper. The 
X, y, z coordinate system is identical with the x, y, z coordi- 
nate systems in Fig. [S] The A and B axes lie in the plane 
spanned by the x and y axes and, in all the structures in 
Figs. [3] (a), (b), and (c), they have the direction of a lattice 
translation. Hence, in all the structures they are perpendicu- 
lar to one another, while the x axis is perpendicular to the y 
axis only in the tetragonal structure. The 2, A, and B axes 
have the same positions in all the three structures in Figs. |3] 
(a), (b), and (c). The point group operations with the indices 
x, y, z, a, and b are related to the depicted x, y, z, A, and B 
axes, respectively. For instance, is a clockwise rotation 
of the lattice through ^ radians about the z axis, and C^a is 
a rotation through vr radians about the A axis, ada = IC2a 
{I is the inversion) is a reflection in the plain containing the 
origin and being perpendicular to the A axis. 



of Imma: One of the four Fe atoms in the unit cell 
Hes at the position pi = ^Ti + + jT^. Applying 
on this atom, e.g., the rotation C'2z, we get an atom at 
the position p2 = — jT^ with the false direction of 
the magnetic moment. Then the associated translation 
t = \T2 + ^Tb, however, moves this atom to the position 
of an Fe atom with the correct direc- 
In the same way, it can 



P3 = •\T2 + iTg 

tion of the magnetic moment 
be shown that the symmetry operation {C22|0ii} leaves 
invariant the directions of the magnetic moments of the 
three other Fe atoms in the unit cell. We get the same re- 
sult for the other two generating elements {C2al50^} and 
{/|^^0}. Consequently, the group Imma is the space 
group of the experimentally determined magnetic struc- 
ture in LaFeAsO. 

In addition, the magnetic structure depicted in Fig. [3] 
(b) is invariant under the anti-unitary operator {if [^^0}, 
where K denotes the operator of time inversion. First, K 
reverses all the magnetic moments (and leaves invariant 
the positions of the atoms) . Then the associated transla- 
tion t — ^Ti + ^T2 moves the Fe atoms to positions with 
the original directions of the magnetic moments. Also 
t leaves invariant the positions of the atoms because it 
is a lattice translation (namely Ts) in the paramagnetic 
lattice depicted in Fig. |3](a). Hence, the magnetic group 
M of the experimentally determined magnetic structure 
in undistorted LaFeAsO as depicted in Fig. [3] (b) may be 
written as 



M = Imma 



{K\^^0}lmma. 



(3) 



Imma, it is sufficient to consider the generating elements 

{C2.|0ii},{C2.|i0i},{/|ii0} (1) 

of Imma as they are given in Table 3.7 of Ref. [l6l |. 
However, in this paper they are written in the coordi- 
nate system defined by Figs. [2] and [3] (b), cf. the notes 
to Table HIl We may easily compare the space group 
Imma with the paramagnetic space group PA/nmm by 
writing the translations of the generating elements ([1]) in 
terms of the basic translations of the Bravais lattice F^ 
of PA/nmm as depicted in Fig. [3] (a). Using the relation 

Ti ^ T1+T3 ] 

T2 ^ -T1+T3 \ Tl ^ Tg (2) 

r3 ^ T2 T3. J 

[as derived from Figs. [3] (a) and (b)] we get the symme- 
try operations {C2z\Jlo}, {Caal^^O}, and {/|001}, all of 
them being in the tetragonal space group PA/nmm, see 
the symmetry operations belonging to point F in Table HI 
Consequently, within the coordinate systems defined in 
Figs. [3] (a) and (b), the group Imma is a subgroup of 
PA/nmm. The important implication is that the posi- 
tions of the atoms of LaFeAsO are invariant under the 
symmetry operations of Imma. 

Now we can show that also the magnetic structure in 
Fig. |3] (b) is invariant under the generating elements ([T]) 



III. 2. The symmetry of the Bloch functions of the 
magnetic band 

The energy band of LaFeAsO denoted in Fig. [1] by the 
bold line is characterized by the representations 

r^, r+; Xi, Xi; M3, Ah; A3, A2; Z+ , ; Ri, R^. (4) 

Folding this energy band into the Brillouin zone of the 
space group Imma of the antiferromagnetic structure in 
the undistorted crystal depicted in Fig. [3](b), the repre- 
sentations (HI) of the Bloch functions transform as 
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see Table IIVI The underlined representations form a 
band hsted in Table El namely band 2 in Tabic |V] (b). 
Hence, the Bloch functions of this band can be unitarily 
transformed into Wannier functions that are 

— as well localized as possible. 
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FIG. 3. Coordinate systems and unit cells of 
three (magnetic) structures in LaFeAsO. For rea- 
sons of clarity, only the Fe atoms and in structure 
(c) also the O atoms are shown, a and c denote 
the lengths of the sides in the tetragonal unit cell. 
The coordinate systems define the symmetry op- 
erations {R\pqr} as used in this paper. They are 
written in the Seitz notation detailed in the text- 
book of Bradley and Cracknell [l^: R stands for 
a point group operation and pqr denotes the sub- 
sequent translation. The point group operation R 
is related to the x,y,z coordinate system as de- 
fined in Fig. [2I and pqr stands for the translation 
t = pTi + qT2 -\- rTs with the basic translations 
Ti being different for the structures (a), (b), and 
(c). The origin of the coordinate systems is fixed 
for all the three structures (a), (b), and (c). 

(a) The paramagnetic structure of LaFeAsO with 
the tetragonal space group PA/nmm (129). 

(b) The antiferromagnetic structure in the undis- 
torted material with the orthorhombic space 
group Imma (74). 

(c) The antiferromagnetic structure in distorted 
LaFeAsO with the "allowed" orthorhombic space 
group Pnn2 (34). This paper proposes the in- 
dicated displacements of the Fe and O atoms in 
z direction which are invariant under the mag- 
netic group M34 = Pnn2 + {K\^^0}Pnn2. Since 
Pnn2 is a subgroup of Imma, the magnetic mo- 
ments may have the same directions as in the 
undistorted structure (b). However, in the mag- 
netic group M34 the moments may also be canted 
within a plain perpendicular to the z axis as de- 
picted in Fig. m 
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— centered at the Fe atoms, 

— and symmetry-adapted to the magnetic group M 
in Eq. ©, 

see the notes to Table |Vl For this reason, this band is 
cahed "magnetic band" related to the magnetic group M . 
The NHM predicts that the electrons of this partly filled 
band may lower their Coulomb correlation energy by ac- 
tivating an exchange mechanism producing a magnetic 
structure with the magnetic group M and with the mag- 
netic moments lying on the Fe sites, i.e., by produ cing th e 
experimentally determined magnetic structure [gI JITIITsI . 
In this sense we say that the symmetry of the Bloch func- 
tions is "compatible" with the experimentally determined 
magnetic structure. 

Table IVl lists all the possible magnetic bands related to 
the magnetic group M in Eq. ([3]). By this table one can 
make sure that the symmetry of the Bloch functions in 
relation ^ corresponds only to band 2 in Table |V] (b) . 
That is, the magnetic moments in LaFeAsO may only be 
situated at the Fe sites. 



IV. STABILITY OF A MAGNETIC 
STRUCTURE 

The magnetic band denoted in Fig. [T]by the bold line 
is related to the magnetic group M in Eq. ([3]). However, 
this structure cannot exist in the undistorted material 
because the space group Imma does not possess suit- 
able representations. This shall be substantiated in the 
following. 

Let be 

M ^ S + {K\t}S (6) 

the magnetic group of a given antiferromagnetic struc- 
ture, with S denoting the space group of this structure 
and jJ^I^} being an anti-unitary operator leaving invari- 
ant the magnetic structure. K still denotes the operator 
of time inversion. 

The operator K reverses the magnetic moments in the 
antiferromagnetic ground state \G), so 

|G) = K\G) (7) 

is the state with the opposite directions of the magnetic 
moments which clearly is different from \G). Within the 
NHM, both states |G) and \G) are eigenstates of a Hamil- 
tonian denoted in Ref. [l7| by H which commutes with 
K. Hence, \G) and \G) belong to an irreducible two- 
dimensional corepresentation D*^ of the group 

M = M + KM (8) 

of the Hamiltonian H, see Ref. 17], Sec. HLC. 

However, when we restrict ourselves to the symmetry 
operations P of the subgroup M of M, then Z?*^ must 



subduce a one- dimensional corepresentation D of M, 
that is, 

P\G) = c • |G) for P e M, (9) 

where |c| = 1. In particular, this Eq. ([9|) is satisfied for 
the anti-unitary operator {Xji}, 

{K\t}\G)^c-\G), (10) 

where still |c| = 1. 

A stable magnetic state \G) may exist if M possesses 
at least one suitable corepresentation D^^ . From this 
condition it follows the 

Theorem IV. 1 A stable magnetic state with the space 
group S can exist if S has at least one one- dimensional 
single-valued representation D 

(i) following case (a) with respect to the magnetic 
group S + {K\t}S and 

(ii) following case (c) with respect to the magnetic 
group S -\- KS. 

The cases (a) and (c) are defined in Eqs. (7.3.45) and 
(7.3.47), respectively, of Ref. [16]. They determine the di- 
mension of the corepresentations of the magnetic groups 
S + {/-iritis' and S -\- KS, respectively, which are de- 
rived 16] from the representation D of S. The one- 
dimensional representation D may satisfy the second con- 
dition (ii) only if D has non-real characters. 

This Theorem II V . 1 1 was proposed in R ef. Il7| and writ- 
ten down in the present form in Ref. [13] • It can be 
understood following the theory of corepresentations as 
given in Sec. 7.3. of the textbook of Bradley and Crack- 
nell [l^: Let D be a representation of S satisfying the 
first condition (i) of the theorem, then we may derive 
from D the one-dimensional corepresentation D*^ of M 
given by Eq. (7.3.45) of Ref. [H (with A = {K\t} and 
N = 1). 

Then we may derive from D^^ the two-dimensional 
corepresentation of M as it is given in Eq. (7.3.17) of 
Ref. [l^, where now A = K. [Bradley and Cracknell as- 
sume that the subgroup G in their Eq. (7.3.11) does not 
contain anti-unitary elements; for the derivation of their 
Eq. (7.3.17), however, these authors do not make use of 
this assumption.] If D satisfies additionally the second 
condition (ii) of the theorem, then the corepresentation 
D^^ has just the required properties, cf. the Auxiliary 
Publication, citation 9, of Ref. [l7| . 

In the present paper the cases (a) and (c) are deter- 
mined (in Tables HIl and IIII[) by the relatively straightfor- 
ward equation (7.3.51) of Ref. (T6j . 

{-\-\S\ in case (a) 
-\S\ in case (b) (11) 
in case (c). 

The sum runs over the symmetry operations B in the 
left cosets C = {K\t}S and C = KS, respectively, of 
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the magnetic groups. x(-B^) is the character of in the 
representation D and \S\ denotes, as usual, the order of 
S. 

Table |lT] shows that the space group Imma does not 
possess one-dimensional representations with non-real 
characters. Consequently, all the one-dimensional rep- 
resentations of Imma follow case (a) with respect to the 
magnetic group Imma + KImma and hence, stable mag- 
netic structures with the space group Imma do not exist. 



V. ALLOWED SPACE GROUPS IN 
ANTIFERROMAGNETIC LaFeAsO 

As shown in the preceding Sec. lIVl the magnetic struc- 
ture depicted in Fig. [H] (b) is not stable because its space 
group Imma does not possess suitable representations. 
However, this structure may be stabilized by a (small) 
distortion of the crystal that turns the space group Imma 
into a new space group possessing representations that al- 
low the formation of a stable magnetic structure. This 
new space group shall be named an "allowed" space 
group in LaFeAsO. In this context we assume that the 
distortions of the crystal are so small that the magnetic 
structure in the distorted material differs hardly from the 
structure depicted in Fig. [3] (b) because the symmetry of 
the Bloch functions of the magnetic band is compatible 
with the magnetic group M in Eq. [31 For that reason we 
assume that the magnetic structure depicted in Fig.[3](b) 
is invariant under an allowed space group, i.e., we assume 
that an allowed space group is a subgroup of Imma. 

In this section we look for all the allowed space groups 
in antiferromagnetic LaFeAsO by considering all the 
orthorhombic primitive and monoclinic primitive space 
groups in Table 3.7 of Ref. [l^ with at least four point 
group elements. This is a laborious task because all the 
conceivable translations and rotations of the coordinate 
systems given in Table 3.7 of Ref. [ll] should be regarded. 

The coordinate system of the magnetic structure in 
LaFeAsO in an orthorhombic primitive or monoclinic 
primitive Bravais lattice is depicted in Fig. [3] (c). The 
positions of the atoms and the magnetic structure are 
invariant under two anti-unitary operations of the form 
{K\t}, namely under {K\^^0} and {K\00^}. Hence, we 
must look for space groups S possessing one-dimensional 
representations following case (c) with respect to the 
magnetic group S + KS and case (a) with respect either 
to S + {K\\\Q}S ov to S + {K\ml}S, see TheoremUvH 

We found the two allowed space groups Pnn2 = VoCl^ 
(34) and PiPiPi = Tp D^ (1 9) a nd sh all present them in 
the following two Sees. ED and ESI In Sec. IvH we will 
show that Pnn2 clearly is the space group of LaFeAsO. 

V.l. The space group Pnn2 

The space group Pnn2 has the orthorhombic-primitive 
Bravais lattice Fq and can be defined by the two gener- 



ating elements 

{CTdalOii} and Kbiiii}, (12) 

cf. Table 3.7 of Ref. [11] and the notes to Table HhI 
The symmetry operations (|12p are given in the coordinate 
system defined in Figs.[2]and[3](c). In order to show that 
Pnn2 is a subgroup of the group Imma (with the Bravais 
lattice Fq) we write the translations of the generating 
elements (|12p in terms of the basic translations of the 
group Imma as given in Fig. [3] (b). Using the relation 

Ti ^ T2+T3 ] 

T2 ^ Ti+Ts \ To ^ ri (13) 
T3 ^ Ti+T2, J 

[as deduced from Figs. [3] (b) and (c)], we get the oper- 
ations {cdall^^} and {(Tdbllll}, both belonging to the 
group Imma, see the symmetry operations belonging to 
point F in Table |lTl Consequently, within the coordinate 
systems of Fig. |3l Pnn2 is a subgroup of Imma. 

In contrast to the group Imma, the group Pnn2 pos- 
sesses at point S one-dimensional representations satis- 
fying the conditions (i) and (ii) of Theorem IIV.II for the 
magnetic groups 

Mu = Pnn2 + {K\^^0}Pnn2 (14) 

and Pnn2 + KPnn2, respectively, see Table IIIIl Fur- 
ther, the little group of S comprises the whole space 
group. Hence, Pnn2 may be the space group of a stable 
magnetic structure. A stable antiferromagnetic ground 
state \G) is basis function of a one-dimensional corep- 
resentation of M34 and \G) and the time inverted state 
\G) = K\G) form basis functions of an irreducible two- 
dimensional corepresentation of the group 

M34 = M34 + KMsi 

which might be determined by Eq. (7.3.17) of Ref. [l^ . 
[As an example, the analogous corepresentation for the 
spin-density- wave state in chromium is explicitly given in 
Ref. [13.] 

V.2. The space group Pi Pi Pi 

A second allowed space group is the group Pi Pi Pi with 
the generating elements 

{CsaliOi} and {CsbjiiO}, (15) 

written again within the coordinate systems of Figs. [2] 
and[3](c). We do not describe this group but only report 
the result that a magnetic structure with the magnetic 
group 

Mi9 = PiPiPi + {X|iiO}PiPiFi (16) 

could be stable in LaFeAsO. However, antiferromagnetic 
LaFeAsO does not possess this space group, see the fol- 
lowing Sec. IVH 
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VI. DISTORTION OF THE CRYSTAL 

The antifcrromagnetic structure in undistorted 
LaFeAsO has the space group Imma. Hence, in the 
undistorted antifcrromagnetic material the nonadiabatic 
Hamihonian [as defined in Eq. (2.15) of Ref. 6]] 
would commute with all the symmetry operations of 
Imma. However, as shown in the preceding sections, 
Imma is not an allowed space group but has two allowed 
subgroups, namely Pnn2 and Pi Pi Pi. Thus, either 
Pnn2 or Pi Pi Pi, but not Imma may be the space 
group of the antifcrromagnetic structure in LaFeAsO. 
Consequently, the nonadiabatic Hamiltonian i/" com- 
mutes with the symmetry operations P of one of the 
allowed space groups S, but does not commute with 
those symmetry operations that belong to Imma, but 
do not belong to S, 



P] = 
P] + 



for 
for 



P G {Imma 



S). 



(17) 



For that reason, the Fe and O atoms in LaFeAsO are 
slightly shifted from their positions in the space group 
Imma. 

In order to show this, we write the eight symmetry 
operations of Imma as given in Table for point F in 
terms of the basic translations Ti, T2, and T3 of the 
orthorhombic primitive lattice Fq in Fig.[3](c), 



{E|000}, 
{C2.|i00}, 

{C2a|0i0}, 
{(Jda\\m, 

{adblOOO}, 



{Em). 



iol}, 



(18) 



{C2. 

{C2a \ 2" 2 
{C2fc|00i}, 

{/|00i}, 

Kioio}, 

{cTdalOii}, 

and {(7db\\\\}, 



cf. relation (|13p . For each symmetry operation from 
Fq we get two operations in F^ because the operation 
in To is a lattice translation in F^ {E is the 

identity) . 

The antifcrromagnetic structure with the space group 
Imma is invariant under the sixteen symmetry opera- 
tions (jl8p now written within the Bravais lattice Fq of 
the structure in Fig. |3] (c). By inspection of Fig. [3] (c) 
we find under these sixteen operations four operations, 
namely 



{i?|000}, 

{CdalOii}, 



{C2.|i00}, 
and {(7db\\\\}, 



(19) 



leaving invariant the distorted crystal as depicted in Fig. [3] 
(c) . These four operations form (together with the trans- 
lations) a group, namely the allowed space group Pnn2. 
Hence, the experimentally observed Q displacements of 
the Fe and O atoms "realize" the space group Pnn2 in the 
sense that they produce the correct commutation prop- 
erties [given by Eq. (IT7|) ] of the nonadiabatic Hamilto- 
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FIG. 4. The two layers of the Fe atoms in the unit cell of 
the distorted antiferromagnetic structure depicted in Fig. [3] 
(c). The z axis is normal to the paper. For clarity, the small 
displacements of the Fe atoms (in z direction) are not indi- 
cated. The X, y, z coordinate system coincides with the x, y, z 
coordinates in Fig. [3] (c). Within distorted LaFeAsO, the 
magnetic moments as depicted in Fig. [3] (c) may be rotated 
by arbitrary angles a and P within a plain perpendicular to 
the z ajds because also this canted structure is invariant un- 
der the magnetic group M34 = Pnn2 + {K\ i iO}Pnn2. a and 
P may be different and both angles are likely small because 
the symmetry of the Bloch functions of the magnetic band 
is compatible with the magnetic group AI of the magnetic 
structure depicted in Fig. Kb). 



nian iJ" and stabilize in this way the magnetic struc- 
ture in LaFeAsO. Together with the anti- unitary oper- 
ator {_fi'|iiO} (which also leaves invariant the distorted 



structure), the symmetry operations (jl9p form the mag- 
netic group M34 in Eq. 

In the magnetic group M34 the magnetic moments may 
be canted by two arbitrary angles a and /3 as depicted 
in Fig. m (while such canted moments are not invariant 
under Imma). Both angles are likely small because the 
symmetry of the Bloch functions of the magnetic band 
is compatible with the magnetic group M [in Eq. ^] of 
the structure depicted in Fig.[3](b). 

For the second allowed space group PiPiPi, on the 
other hand, we have not found any simple distortion of 
the lattice realizing this group. 



VII. CONCLUSIONS 

Within the nonadiabatic Heisenberg model (NHM), 
the situation in the magnetic band of LaFeAsO (denoted 
in Fig. [1] by the bold line) is characterized by two dif- 
ferent group-theoretical phenomena. On the one hand, 
the magnetic band (i.e., the symmetry of the Bloch func- 



8 



tions of this band) is related to the magnetic group M 
of the experimentally observed magnetic structure 
in LaFeAsO and, on the other hand, any structure with 
the space group M is unstable. These two phenomena 
are discussed in separated subsections IVII.ll and IVII.2I 



VII. 1. Nonadiabatic condensation energy 

The nonadiabatic Heisenberg model (NHM) is defined 
by three postulates on the Coulomb correlation energy in 
narrow, partly filled bands (61]. A direct consequence of 
these postulates is the existence of the "nonadiabatic con- 
densation energy" l^E defined by Eq. (2.20) of Ref. [i]. 
The electrons at the Fermi level may lower their Coulomb 
correlation energy by Ai? by occupying an atomic-like 
state as defined by Mott [2^ and Hubbard |2ll] : the elec- 
trons occupy localized states as long as possible and per- 
form their band motion by hopping from one atom to 
another. In the present approach, the atomic-like state 
is consistently described in terms of symmetry-adapted 
and optimally localized Wannier functions and is pre- 
cisely defined by an equation in the nonadiabatic system, 
namely by Eq. (2.19) of Ref. Q. In case the considered 
partly filled band is a magnetic band related to a mag- 
netic group M (as considered in this paper), the electrons 
of this band may gain the energy /S.E by activating an ex- 
change mechanism producing a magnetic structure with 
the magnetic group M [13, [3 • 

In former papers we could show that evidently the 
magnetic states in Cr [l3j, Fe La2Cu04 fl^, and 
YBa2Cu306 are connected with narrow, partly filled 
magnetic bands in the band structures of the respective 
materials. That is, there is evidence that the nonadia- 
batic condensation energy /S.E is responsible for the oc- 
currence of the magnetic states in these materials. 

Also LaFeAsO possesses a magnetic band related to 
the magnetic group M = Imma + {K\^^0} of the ex- 
perimentally observed magnetic structure. This finding 
suggests that, first, the electrons of the magnetic band 
of LaFeAsO perform the atomic-like motion as defined 
within the NHM by Eq. (2.19) of Ref. @ and that, sec- 
ondly, the nonadiabatic condensation energy AE is re- 
sponsible for the magnetic state. 



VII. 2. Stable magnetic structures 

According to Theorem lIV.il not all the magnetic struc- 
tures may be stable within the NHM. This theorem led al- 
ready to an understanding of the really existing magnetic 



groups in Cr [13, La2Cu04 [H, and YBaaCugOe ^. In 
the present paper, this theorem is the key for an un- 
derstanding of the distortion (in addition to the magne- 
tostriction) of antiferromagnetic LaFeAsO. 

Theorem IIV.II defines "allowed" space groups belong- 
ing to stable magnetic structures. Though the mag- 
netic band of LaFeAsO is related to the magnetic group 
M = Imma + {K\^^0}, the space group Imma is not 
allowed. On the other hand, allowed space groups in 
LaFeAsO are the two groups Pnn2 and Pi Pi Pi. Con- 
sequently, the magnetic structure may be stabilized by 
a (slight) change of the atomic sites in such a way that 
the space group Imma is turned into one of the allowed 
space groups. Briefly speaking, the structural distortion 
"realizes" the space group allowed in LaFeAsO. 

The displacements of the Fe and O atoms as proposed 
in this paper are depicted in Fig. |3] (c). They have the 
space group Pnn2 and realize exactly this allowed space 
group. These displacements coincide in essence with 
the theoretically determined (2^ and experimentally ob- 
served [3| structural distortions of LaFeAsO. 

However, we could not confirm the space groups 
Cmma and P112n as they have been reported for the 
space CTOup of antiferromagnetic LaFeAsO in the undis- 
torted [l| and the distorted Q crystal, respectively. Nev- 
ertheless, using the space group P112n, Clarina de la 
Cruz et al. (3| discovered the important displacements of 
the Fe and O atoms in z direction. The directions of the 
displacements depicted in Fig. [3] (c) do not essentially 
differ from those given in Table 2 of Ref. 3 : while they 
coincide absolutely within the Fe or O layers. Table 2 
of Ref. suggests that the displacements are periodic 
with T3/2. As depicted in Fig. [3] (c), the displacements 
proposed in this paper alternate their directions at an 
interval of T3/2. However, in the space group Pnn2 the 
amounts of the displacements in up and down direction 
may be principally different because they are not con- 
nected by symmetry operations. 
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Appendix: Tables 



TABLE I. Character tables of the single-valued irreducible representations of the space group PA/nmm = T^Dlh (129) 
tetragonal paramagnetic LaFeAsO. 

r(ooo), z(ooi) 

{Q-JOiO} {C2j;|OiO} {C25IOOO} {S'4+IOiO} {aJOiO} {tTdblOOO} 

{ig|ooo} {ctAhm {c^Ahm {C2a\khQ} Ulooo} KlUo} {^rJloo} KI|oo} Wall^o} 



1 1 


^1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


r+ 




1 


1 


1 


-1 


-1 


1 


1 


1 


-1 


-1 


r+ 
^ 3 ' 


^3+ 


1 


1 


-1 


1 


-1 


1 


1 


-1 


1 


-1 




zt 


1 


1 


-1 


-1 


1 


1 


1 


-1 


-1 


1 




zt 


2 


-2 











2 


-2 













z^ 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


r2-, 


Z2 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


1 


1 




Zz 


1 


1 


-1 


1 


-1 


-1 


-1 


1 


-1 


1 


r4-, 


z^ 


1 


1 


-1 


-1 


1 


-1 


-1 


1 


1 


-1 




zi 


2 


-2 











-2 


2 












M(llO) 



{<Trf„|iiO} {a,,„|ifO} {/|oin} 
{E\()m) {L'lOiO} {C2.IH0} {':"'2.|jrt"} I'^dblOOO} |ad6l010} {iiooo} 



Ml 


2 


-2 


2 


-2 


2 


-2 





M2 


2 


-2 


-2 


2 











M3 


2 


-2 


2 


-2 


-2 


2 





M4 


2 


-2 


-2 


2 












M(iiO) (continued) 

Kiofo} {c+|iio} {C2,|o|o} {54- I i 10} 

KlOiO} {Q-jOfO} {C2,|0i0} {S+|OiO} 
{a,|iiO} {6-26 1000} {Cat, iOlO} K|ilO} {Q-JOiO} {C2.||00} {S4+|0|0} 
{C2a\hl0} {C2a\kk0} KI^OO} {^4+1^00} {^2.1^10} {^4'J^oo} 



Ml 























M2 





2 


-2 














Ms 























M4 





-2 


2 
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^(555) 







{o-da ifiO} {o-daliil} {/|001} 




{£1000} {£1001} {C2.|iiO} {C2.|iil} {ddblOOO} {ad^lOOl} {/|000} 




2 -2 


2-22-20 


A2 


2 -2 


-2 2 


A3 


2 -2 


2 -2-2 2 




2 -2 


-2 2 






yl(ili) [continued) 






{cjy\Q\l} {Q+JiOl} {C2,|0il} {SrJiOl} 






{a,|0i0} {Q-JOil} {C72,|0i0} {S4+|OiO} 


K|iiO} {C26IOOO} 


{CatlOOl} {a,|i01} {C4.|OiO} {C2.|i00} {S+JOil} 


Kjiil} {C2a|Ul} {C2a\k\Q} W|iOO} {Q+jiOO} {C2.[i01} {^rJiOO} 


Al 










2 


-2 


A3 








A4 


-2 


2 



ii(Oii) 

{C2,10ll} {C72.[|00} {C2.|ii0} {'^Ahk^} {/|000} {a,|OiO} 
{£|000} {£|001} {C2y|0i0} {C2.|^01} {C2.|^^l} {t^^l^OO} {a.|^01} {a-Ahk^} {^|ooi} Klo|i} 
i?i 2 -2 2 -2 
R2 2 -2 -2 2 

X(OiO) 

{C2y|0|0} {C2.|i00} {C2.|ii0} KllfO} {^1000} {(7,,|OiO} 

{£|000} {£|010} {Ca^lO^O} {C2x||l0} {C2z|Mo} Kl^oo} Kl^io} KlHo} {-^1010} Klofo} 
Xi 2 -2 2 -2 

X2 2 -2 -2 2 



Notes to Table U 

(i) The space group operations are related to the coordinate systems in Figs. [2] and [3] (a). 

(ii) The character table is determined from Table 5.7 of Ref. The origin of the coordinate system used in Ref. ;lf| is 
translated into the origin used in this paper by to = —\Ti — \T2. That is, the symmetry operations Phc given in Table 
5.7 of Ref. [1^ are changed into the operations Pp used in this paper by the equation 

Pp = {E\tQ}Pi,{E\ ~ io}, 



where E is the identity operation, cf. Eq. (3.5.11) of Ref. [IGj (which is related to the opposite translation —to). 



12 



TABLE II. Character tables of the single-valued irreducible representations of the orthorhombic space group Imma = r]l,D2h 
(74) of the antiferromagnetic structure depicted in Fig. [3](b). 





K 


{£|000} {C2.|0ii} {C2a|| 


r(ooo) 

Qi} {C2i,|iio} {/IMo} K 




IO55} Wblooo} 




(a) 


1 


1 


1 


1 


1 


1 


1 


1 




(a) 


1 


-1 


1 


-1 


1 


-1 


1 


-1 




(a) 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


r+ 


(a) 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


rr 


(a) 


1 


1 


1 


1 

i 


1 

-i 


-1 


-i 


-i 




(a) 


1 


-1 


1 


-1 


-1 


1 


-1 


1 




(a) 


1 


1 


-1 


-1 


-1 


-1 


1 


1 




(a) 


1 


-1 


-1 


1 

i 


-1 


1 




-i 




K 


{£1000} {C2.|0ii} {C2a\\ 


X (222) 

oi} {C26|iil} {/ 


1 1 1 T r 


l|oi} W 




o-dblOOO} 




(a) 


1 


1 


1 


1 


1 


1 


1 


1 




(a) 


1 


-1 


1 


-1 


1 


-1 


1 


-1 




(a) 


1 


1 


-1 


-1 


1 


1 


-1 


-1 


xt 


(a) 


1 


-1 


-1 


1 


1 


-1 


-1 


1 


x^ 


(a) 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


X2 


(a) 


1 


-1 


1 


-1 


-1 


1 


-1 


1 


x^ 


(a) 


1 


1 


-1 


-1 


-1 


-1 


1 


1 


XI 


(a) 


1 


-1 


-1 


1 


-1 


1 


1 


-1 



i?(iOO) S{^0^) 





K 


{£1000} {C2a 


I^Oi} {J| 


i-|0} {Oda 


loii} 




K 


{£1000} {C2b 


1^0} {/| 


Ml} {'^dilOOl} 


Rt 


(a) 


1 


1 


1 


1 


St 


(a) 


1 


1 


1 1 


Ri 


(a) 


1 


1 


-1 


-1 




(a) 


1 


1 


-1 -1 


Rt 


(a) 


1 


-1 


1 


-1 


St 


(a) 


1 


-1 


1 -1 


R2 


(a) 


1 


-1 


-1 


1 


S2 


(a) 


1 


-1 


-1 1 



Notes to Table HI] 

(i) The points T and W are not listed because they possess only one two-dimensional representation T\ and W\ , respectively. 

(ii) The character tables are determined from Table 5.7 in Ref. The origin of the coordinate system used in Ref. [3| 
is translated into the origin used in this paper by to = ^Ta. That is, the symmetry operations Pbc given in Table 5.7 of 
Ref. are changed into the operations Pp used in this paper by the equation 

Pp = {E\to}Pt,{E\ - to}, 

where E is the identity operation, cf. Eq. (3.5.11) of Ref. [iQ] (which is related to the opposite translation —to). 

(iii) The space group operations are related to the coordinate systems in Figs. [2] and |3](b). The x, y, and z axes have the 
same orientations as in the tetragonal structure in Fig. [3] (a) . In this way, Imma becomes a subgroup of the tetragonal 
group P4:/nmm. As a consequence, the xyz coordinate system used in Ref. [l^ is rotated anti-clockwise with respect 
to the basic translations through j radians about the z axis. Thus, the x and y axes in Ref. [l^ become the B and 
A axes, respectively, in this paper, see Fig. [S] Consequently, the point group operations belonging to the space group 
operations calculated by the above equation are renamed in a second step. For instance, the operation C2x in Table 
5.7 of Ref. [Igl turns into the operation C26. 

(iv) K stands for the operator of time-inversion. The entries below K specify whether the related representation follows, with 
respect to the magnetic group Imma + KImma, case (a), case (b), or case (c) when they are given by Eqs. (7.3.45-47) 
of Ref. [li|. 

(v) The cases (a), (6), and (c) are determined by Eq. (|lip . 

(vi) All the one-dimensional representations of the space group Imma are real and, hence, follow case (a) with respect to 
Imma + KImma. Therefore, stable magnetic structures with this space group do not exist, see Theorem lIV.il 
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TABLE III. Character tables of the single-valued irreducible representations of the orthorhombic space group Pnn2 = roC22 
(34) of the antiferromagnetic structure depicted in Fig. [3](c). 

r(ooo) 



K 


{/CliiO} {i^lOOi} {£|000} {C2.|i00} {aaa 


|0| 


5} {o-di>llM} 


Ti (a) 


(a) (a) 1 1 


1 


1 


(a) 


(a) (a) 1 1 


-1 


-1 


r2 (a) 


(a) (a) 1 -1 


1 


-1 


r4 (a) 


(a) (a) 1 -1 


-1 


1 



Z(OOi) 

K {K\^^0} {K\00^} {g|000} {ada\Ok^} {E\001} Wa|OM} {^2.1^00} WtlH^} {Ca.l^Ol} Wdt\kk'^} 





(c) 


(c) 


(c) 


1 


i 


-1 


-i 


1 


i 


-1 


-i 


Z2 


(c) 


(c) 


(c) 


1 


-i 


-1 


i 


1 


-i 


-1 


i 


Z-i 


(c) 


(c) 


(c) 


1 


i 


-1 


-i 


-1 


-i 


1 


i 


Za 


(c) 


(c) 


(c) 


1 


-i 


-1 


i 


-1 


i 


1 


-i 



S(iiO) 



K 


{K\U 


0} {K\00^} {£1000} {ado 


jQii} {£1010} {aaa 


|0|i} {C2.|il0} {aat 


|i 1 1} {C2.|i00} {a,,ii|i} 


Si (c) 


(a) 


(c) 


1 


i -1 


-i 1 


i -1 -i 


S2 (c) 


(a) 


(c) 


1 


-i -1 


i 1 


-i -1 i 


S3 (c) 


(a) 


(c) 


1 


i -1 


-i -1 


-i 1 i 


S4 (c) 


(a) 


(c) 


1 


-i -1 


i -1 


i 1 -i 



Notes to Table Iml 

(i) In addition to F only the points Z and S are listed because only Z and 5* possess one-dimensional non-real representations. 

(ii) The character tables are determined from Table 5.7 in Ref. [l^. The origin of the coordinate system used in Ref. [lH| 
is translated into the origin used in this paper by to = jT2. That is, the symmetry operations Pbc given in Table 5.7 of 
Ref. [l6|] are changed into the operations Pp used in this paper by the equation 

Pp = {E\to}Pic{E\ - to}, 

where E is the identity operation, cf. Eq. (3.5.11) of Ref. [161 ) (which is related to the opposite translation —to). 

(iii) The space group operations are related to the coordinate systems in Figs. [2] and [3] (c). The x, y, and z axes have the 
same orientations as in the tetragonal structure in Fig. |3] (a). In this way, Pnn2 becomes a subgroup of the tetragonal 
group P4:/nmm. As a consequence, the xyz coordinate system used in Ref. [l^ is rotated anti-clockwise with respect 
to the basic translations through j radians about the z axis. Thus, the x and y axes in Ref. [l^ become the B and 
A axes, respectively, in this paper, see Fig. (2] Consequently, the point group operations belonging to the space group 
operations calculated by the above equation are renamed in a second step. For instance, the operation in Table 
5.7 of Ref. [la] turns into the operation adt- 

(iv) K stands for the operator of time inversion. The entries below K and the operators {K\t} specify whether the related 
representation follows, with respect to the magnetic groups Pnn2 + KPnn2 and Pnn2 + {K\t}Pnn2, respectively, case 
(a), case (6), or case (c) when they are given by Eqs. (7.3.45-47) of Ref. fl^. 

(v) The cases (a), (6), and (c) are determined by Eq. (|lip . 

(vi) The little group of point S comprises the whole space group Pnn2 and S possesses one-dimensional representations 
following case (c) and case (a) with respect to Pnn2 + KPnn2 and Pnn2 + {K\j^Q}Pnn2, respectively. Consequently, 
a magnetic structure with the space group Pnn2 can be stable, see Theorem lIV.il 
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TABLE IV. Compatibility relations between the Brillouin zone for tetragonal paramagnetic LaFeAsO and the Brillouin zone 
for the orthorhombic antiferromagnetic structure depicted in Fig. [3] (b). 

r X 



r+ r+ r+ r+ r+ r- r- rj r- Kh Ah M3 M4 

^3 r+ + r+ r- r3 r- r- r,- + rj x+ + x^ x+ + x^ x+ + x^ x+ + x^ 



X 



Z+ Z+ Z+ Z+ Z+ Z- Z^ Z^ Zl Z- A-v M M A4 

xt x+ x+ x+ x+ + x+ X- X.- X- X- x^ + X4- r+ + r- r+ + r,- + r- r+ + rj 



T 



Ri R2 Xi X2 

Ti Ti Ti Ti 



Notes to Table HV] 

(i) The antiferromagnetic structure in Fig. [3] (b) has the space group Imma. 

(ii) The Brillouin zone for the orthorhombic space group Imma lies within the Brillouin zone for the tetragonal space group 
PA/nmm. The points F and A in the tetragonal Brillouin zone are equivalent to the point F in the orthorhombic Brillouin 
zone; M and Z are equivalent to X, and R and X are equivalent to T. 

(iii) The upper rows list the representations of the little groups of the points of symmetry in the Brillouin zone for the 
tetragonal paramagnetic phase. The lower rows list representations of the little groups of the related points of symmetry 
in the Brillouin zone for the antiferromagnetic structure. 

The representations in the same column are compatible in the following sense: Bloch functions that are basis functions 
of a representation Di in the upper row can be unitarily transformed into the basis functions of the representation given 
below Di. 

(iv) The compatibility relations are determined in the way described in great detail in Ref. [2^ . 

(v) The representations are labeled as given in Tables U and [ITl respectively. 
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TABLE V. Single-valued representations of all the magnetic bands related to the space group Imma = T^D^t (74) with the 
Wannier functions being centered (a) at the La atoms, (b) at the Fe atoms, (c) at the As atoms, and (d) at the O atoms. 

(a) La 

T X R S T W 

Band 1 2r+ + 2X+ + 2X~ 2R+ + 2R~ 2S~ + 2S} 2Ti 2Wi 
Band 2 2rJ + 2rj 2X^ + 2X1 '^Rf + 2i?2 "^S^ + 23^ 2Ti 2Wi 
Band 3 2r+ + 2rj" 2X^ + 2X^ 2R^ + 2i?J 25^ + 2S2 2Ti 2Wi 
Band 4 2rj + 2r2 2X/ + 2X2 2R^ + 2Rt 2gf + 2S} 2Ti 2Wi 



(b) Fe 

T X R S T W 

Band 1 r+ + r+ + r- + X+ + X+ + X- + X,J R+ + R^ + R+ + R^ S+ + S- + S+ + 2Ti 2Wi 
Band 2 r+ + r+ + F^ + X+ + X+ + X' + X^ R+ + R- + Rj + R^ S+ + + S+ + g~ 2Ti 2Wi 



(c) As 

F X R S T W 

Band 1 2F^ + 2F3 2X+ + 2X^ 2R+ + 2i?2 ^S' + 2St 2Ti 2Wi 
Band 2 2Fj + 2FJ 2X^ + 2X4" 2R+ + 2i?2 25'^ + 25'2" 2Ti 2Wi 
Band 3 2F+ + 2Fj; 2X+ + 2X[' 2R^ + 2R+ 2S^ + 23^ 2Ti 2Wi 
Band 4 2r+ + 2r2 2X+ + 2X2" 2R:[ + 2R+ 2S^ + 2g+ 2Ti 2Wi 



(d) O 

F X R S T W 

Band 1 F+ + F+ + Fj: + F2" A:+ + X+ + X' + X^ Rt + R^ + Rt + R2 + ST + 5*2+ + 2Ti 2Wi 
Band 2 F+ + F+ + F3 + F4 X+ + X+ + X^ + X^ R+ + R{ + R} + R^ S+ + + S+ + ^2' STi 21^1 



Notes to Table |V] 

(i) The antiferromagnetic structure of LaFeAsO depicted in Fig. |3] (b) has the space group Imma and the magnetic group 
M = Imma + {K\^^0}lmma with K denoting the operator of time-inversion. 

(ii) Each row defines one band consisting of four branches, because in each case there are four atoms in the unit cell. 

(iii) The representations are given in Table [III 

(iv) The bands are determined by Eq. (23) of Ref. 

(v) Assume a band of the symmetry in any row of theses Tables FVl (a), (b), (c), or (d) to exist in the band structure of 
LaFeAsO. Then the Bloch functions of this band can be unitarily transformed into Wannier functions that are 

— as well localized as possible; 

— centered at the assigned (La, Fe, As, or O) atoms; 

— and symmetry-adapted to the space group Imma. 

(vi) Eq. (23) of Ref. [3 makes sure that the Wannier function may be chosen to be symmetry-adapted to the space group 
Imma. In addition, there exists a Matrix N satisfying both Eqs. (26) (with {7^11^0}) and (32) of Ref. 19] for all 
the bands listed in this table. Hence, the Wannier functions may be chosen symmetry adapted to the magnetic group 
M = Imma + {K\\\Q}Imma. 



